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Abstract 

It is well-known that the transport properties of single-electron transistors with a superconducting 
island and normal-conducting leads (NSN SET) may depend on whether or not there is a single 
quasiparticle on the island. This parity effect has pronounced consequences for the linear transport 
properties. Here we analyze the thermopower of NSN SET with and without parity effect, for 
entirely realistic values of device parameters. Besides a marked dependence of the thermopower on 
the superconducting gap A we observe an enhancement in the parity regime which is accompanied 
by a dramatic increase of the thermoelectric figure of merit ZT. The latter can be explained within 
a simple re-interpretation of ZT in terms of averages and variances of transport energies. 

PACS numbers: 72.15.Jf, 73.23.Hk, 74.45. +c 



1 



I. INTRODUCTION 



The properties of transport through small conducting islands have been investigated ex- 
tensively during the past years. Electric current in such devices flows due to tunneling of 
single electrons and is subject to the so-called Coulomb blockade effect |l| which is charac- 
terized by a new energy scale, the capacitive charging energy Eq of the island (see below). 
In the recent past, substantial attention has been devoted also to thermoelectric effects in 



single-electron devices j2|-|l7|. Single-electron devices are interesting candidates for thermo- 



electric applications as it has long been known that dimensiona 
dynamics may lead to an enhanced thermoelectric efficiency 18 



reduction of the electron 



19|. 



While an immense amount of work has been done to investigate thermopower for quantum 
dots, surprisingly little is known about the thermoelectric effects in single-electron transistors 
(SET) with superconducting electrodes. In particular, SET with superconducting islands 
are interesting as they may exhibit the parity effect where a single unpaired quasiparticle 



determines the macroscopic thermodynamic properties of the island electrode 
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2lj as 



well as the cur rent- volt age characteristics of SET with normal-conducting electrodes and 



a superconducting island (NSN SET) 22l424| . The parity effect can be observed below 



a crossover temperature T* ^ A/8 for typical system parameters (here A is the energy 
gap of the superconductor). One may expect that the peculiar combination of properties 
like several competing energy scales {Ec, A, T*) and the presence of a singularity in the 
quasiparticle spectrum gives rise to interesting behavior in the thermoelectric response of 
such systems. 

In Ref. jlll the thermopower for an NSN SET was studied for A < Ec and T > T*. 
Even for this regime without parity effect, interesting oscillations of the thermopower as 
a function of the electrostatic island potential and their strong dependence on the ratio 
A/Ec were predicted. In this article, we investigate the thermopower of NSN SET for 
temperatures below the crossover temperature. We find that the interplay of energy scales, 
Coulomb-blockade and parity effects, and the peculiarity of the electronic spectrum lead 
to a rich variety of features in the thermopower S. Most intriguingly, however, for certain 
gate voltages this system displays a dramatic enhancement of the thermoelectric efficiency 
quantified by the figure of merit ZT. 

This paper is organized as follows. First we introduce the setup and the theoretical 



methods that are used to describe transport in such systems. Subsequently we briefly review 
the parity effect in SET with superconducting islands. We then turn to describe the results 
for the thermopower, and to interpret them in terms of average transport energies. Finally 
we discuss the surprising results for the flgure of merit ZT which we explain in the frame of 
a simple re-interpretation of this quantity. 



II. THE TRANSISTOR SETUP AND MASTER EQUATION 

In an NSN SET, a superconducting island with a small electrostatic capacitance C is 
connected via tunnel junctions to two normal-conducting leads (cf. Fig. 1). The correspond- 
ing charging energy Ec = e^/{2C) ^ ksT is large compared to the temperatures under 
consideration (here, e > denotes the elementary charge). The conductances of the tunnel 
junctions are assumed to be small compared to e^//i, so sequential tunneling dominates and 
cotunneling effects may be neglected. 
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FIG. 1. The NSN SET consists of a superconducting island (S) which is coupled to two normal- 
conducting leads (N) via tunnel barriers. The electrostatic potential of the island can be controlled 
by the gate voltage Vg. There may flow a current through the system due to the bias voltage V or 
a temperature difference AT = Ti —T^ between the two leads. In order to detect the thermopower 
S = —V/AT as a function of the gate voltage Vg the bias V is adjusted such that the corresponding 
current exactly cancels the current which arises due to the temperature difference. 
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The electrostatic potential of the island can be controlled by means of an external po- 
tential Tlx oc Vg due to the gate voltage Vg, and the electrostatic energy of the setup with n 
excess electrons on the island can be expressed as 

Enin^) = Ec {n^ - 2nn^) . (1) 

The energy cost for adding a single electron to the island is Un{nx) = En+i{nx) — En{nx) 
while to add two electrons, an energy En+2{f^x) ~ En{nx) = Un+i{nx) + Uni^x) is required. 
A current may flow in the device if a bias voltage AV^ = Vr — Vi and/or a temperature 
difference AT = Ti — is applied. Throughout this work we consider the linear-response 
regime, that is, \AV/Ec\ < 1 and |Ar/T| < 1. 

Linear transport of charge and heat is conveniently described in terms of the equations 



I, 




(2) 



which relate the charge and heat current response, Jg and Ig, to an applied potential and 
temperature difference. Here, Gy is the linear (charge) conductance, and the thermal con- 
ductance K is defined via Iq = kAT for Je = 0, i.e., k, = K — GyTS"^. The relevant quantity 
for the thermoelectric response is the thermopower S, given by S* = —AV/AT = Gt/Gv- 

In order to calculate the the conductances Gy, Gt and the thermopower S we employ a 
master-equation formalism. According to the orthodox theory l| charge and heat current 
through the system can be written as 

/e = -e 5^ 5^ P„ [vr''-' - r;^"+^] (3) 

n i=l,2 

= E E [?"^"~^' - g"^"^^] (4) 

n i=l,2 

where P„ is the stationary probability for finding n excess electrons on the island, F"^""-^ 
is the tunneling rate of j electrons from the island to the right lead, and T"'^""'"-' denotes 
the tunneling rate of j electrons from the right lead to the island. Correspondingly, g"^""-' 
is the energy transfer rate in j-electron tunneling to the right lead whereas denotes 
the energy transfer rate in a j-electron tunneling event from the right lead to the island. We 
consider only sequential tunneling and neglect co-tunneling events. 

In Eqs. ([3]),(jl]) we have taken into account the possibility of single-electron tunneling 
{j = 1) and coherent two-electron tunneling (j = 2) ^]. The rates for the latter process 



are given by 

rn^n±2( (±2)^ = ^ ^ , ^ = / ,r (5) 

('+21 

where Ga,/? G^yi,r are the Andreev conductances in the left and the right junction and , 
€r are the energies which are dissipated in a two-electron transfer, e.g.^ = ujynx) + 
M„+i(nx) — 2eV5. The energy transfer rate for two-electron tunneling is then obtained as 

^„^„±2^^(±2)^ = ef^^rr-^^lef^^) (6) 

where the reference point is the Fermi level of the leads. 

The rates for single-electron transitions are given by a sum of two contributions. On 
the one hand, we have the standard expressions for tunneling, e.g.^ from a normal to a 
superconductor |26l | 

r,(e) = ^2 ^E-^===[h{E-e)U-E) + U{-E-e)UE)] .^ = l.r (7) 

(where fi{x) = 1/(1 + exp (x/Ti)) denotes the Fermi function with the appropriate temper- 
ature Ti, Tr, or T for the leads or the island, respectively). On the other hand, t 



lere is the 



escape rate of a single unpaired quasiparticle whose energy equals that of A 
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with the normal-electron density of states per spin on the island z/isi. Note that there is 
also a corresponding recombination rate. In many cases these escape rates can be neglected, 
however, at very low temperatures they may exceed the subgap tunneling rate originating 
from thermally excited quasiparticles in the superconductor. In that case, they produces a 
different macroscopic behavior of the system, depending on whether the total charge number 
on the island is even or odd. Thus we have 

pn->„±i. (±l)^ _ /r.(er^) forneven 
\Ti[el ') for n odd 

where, e.g., e[^^^ = Un{nx) — eVi. Finally, the corresponding energy transfer rates for single- 
electron tunneling are found from 

n(6) = -^2j^ dE-^====[iE-e)ME-e)U-E) + {-E-e)M-E-e)UE)] 



and 



(where again the reference point is the lead Fermi level) such that we have 

^ l^li^'') forneven 

\m^'^)+ne^'^) fornodd ^ ^ 

In order to calculate the currents ([3]), (jlj) we solve the stationary master equation for the 
probabilities Pn 

^ = = J^p.r'^-^-p.r"-^ 

for an applied (small) bias voltage AV^ or a (small) temperature difference AT. If a voltage 
or temperature difference is applied the island is, strictly speaking, in a non-equilibrium 
state. For the given physical situation it is reasonable to neglect the non-equilibrium part 
of the distribution function (which we have already done by writing the transfer rates in 
the form above). For the temperature T of the island we assume the arithmetic mean 
T = (T; + Tr)/2. We note that there is an independent test of the numerical calculation by 
checking the Onsager relation Gt = M/T for the coefficients in Eq. ([21) which is obeyed to 
a high accuracy by our method. 

III. THE PARITY EFFECT 

As mentioned above, the parity effect arises as a macroscopic manifestation of the parity 
of the electron number in a superconductor, i.e., different behavior dependin g on whether 



the total electron number is even or odd. The effect was predicted in Ref. 20j and first 



observed by Tinkham et. al. 



2l|. 



In an even-number superconductor at T = all electrons near the Fermi level are bound in 
Cooper pairs and there is not a single unpaired quasiparticle left. This system has an energy 
gap 2A. On the other hand, adding one electron results in a single quasiparticle excitation 
which which does not have an excitation gap. However, the response to an external field 
is drastically reduced as it is caused just by a single electron. It is intuitively clear that 
the difference between the two parities will fade away as soon as there are more thermal 
quasiparticles. This defines the criterion for a crossover temperature T* beyond which even- 
odd differences disappear for an isolated superconductor such as the island in the NSN SET: 

6 



It is the temperature at which there is on average one thermally excite quasiparticle and it 
is defined by 

T* = (11) 

where Ncs{T) = z/isiv^27rTA can be viewed as an effective number of accessible quasiparticle 
states at the temperature T. Superconducting islands in SET are made of aluminum, and 
for typical parameters (cf. Refs. 2ll-l23|) one has A^ofr ~ 10^ and T* ~ 250mK. 



Keeping in mind that an SET has the characteristic energy scale Ec, there arise four 
interesting transport regimes: For the gap and the charging energy we may have A < Ec 
or A > Eq- These two cases may be studied for T > T* (no parity effects) or in the parity 
regime T < T*. Interestingly, with our numerical method we have a choice for studying 
the system with parity effects. For computations at high temperatures T > T* it does 
not matter whether or not the two-electron and escape rates are included, they do not 
give any observable effect. On the other hand, for T < T* even-odd differences cannot be 
observed without including these rates. Therefore, including or not including these rates 
in calculations below T* helps us to identify the contribution due to the 'parity-generating' 
processes. 



IV. THERMOPOWER OF AN NSN SETUP 



Before we discuss our results we briefly recall an idea due to Matveev 27| (cf. also |ll|) 
to interpret the thermopower as an average energy (^) at which current is transported in the 
voltage-biased system (under linear transport conditions). This is an intuitive and powerful 
method which will use throughout our work to interpret our results. 

The idea can be understood by taking into account that the current Jg through a device 
can be written as I^, = —e /[/i(0 ~ /r(0]^(0'^^ where fi{C,) and /r(0 denote the Fermi 
functions in the left and right lead for the respective potential and temperature (the energies 
^ are taken with respect to the chemical potentials). Here w(^) includes all other quantities 
such as density of states in the leads and transparency of the device at the energy C,- By 
noting that Gt = dlq/dTi and Gy = dl^/dV one flnds 



/(-|)u>(f)df 



and consequently for the thermopower 

S = . (12) 

That is, up to a factor the thermopower measures directly this average energy (,^) . 

Let us turn now to our calculations for the thermopower of NSN SET. For the numerics 
we measure energy and inverse time in units of £'c, and conductance in units of e^. Then, the 
free parameters are the ratios of the conductances G^, G,., Ga,i^ Ga.t^ and the escape rates 
7i, The realistic parameter values which enter our calculations are: Gi = = (50kf2)~^ 
Ga,i = Gji^r = 5 ■ 10^^i7^^, Eq = 100/ieV, 7; = 7r = lO^s^^. The experiments in Ref. 
have been carried out with device parameters close to these values. 
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A. A<Ec 

In Fig. 2 we show the results for calculations of the thermopower with and without parity 
effects. The functional dependence S{nx) in the latter case has been explained analytically 



in Ref. fll|. As we have T <^ A, Ec it is sufficient to include two (or at most three) charge 
states in the considerations in order to understand the behavior of S{nx). Here we focus 
on a brief discussion of the additional features in the parity regime T < T*. The key to 
understand the functional dependence of S{nx) without even-odd effects is that at the points 
Un = 0, that is, when takes half-integer values tunneling to the island n n + 1 and 
from the island ri + 1 — )■ n occurs with equal probability. By using simple arguments for the 



probabilities P„ and Pn+i the following equation has been derived in Ref. 



ill 



S = --^ ( Uo{n^) - A tanh 



2T 



(13) 



which governs the behavior of S{nx) in the intervals from Un{nx) = to the zeros of S{nx) 
which are closest to those values. 

By inspecting Fig. 2 we note that the essential difference introduced in the parity regime 
is that the purple curve is shifted from the half-integer values of Ux towards the closest 
even number. As the slope of the curve does not change in this shift, the consequence is 
that the maximum absolute value of S increases, compared to the case without parity. It 
is not difficult to describe this behavior analytically by repeating the arguments that lead 
to Eq. (fT3|) . now taking into account that the recombination rate of quasiparticles on the 

8 



A=0.5E^, T=0.03E^ 




FIG. 2. Thermopower S{n^) of the NSN SET with A = O.bEc and temperature T = OmEc- The 
green dashed hne represents the thermopower without two-electron tunneling and escape rate and 
therefore displays no even-odd differences. The purple solid line displays the result including all 
processes. It shows clear 2e periodicity and an enhancement of the thermopower for certain gate 
charges. 



island is no longer oc exp(— A/T). In this range of there is an unpaired quasiparticle 
on the island whose recombination rate provides the dominating contribution to the charge 
current. This rate is oc (e~^/^ + t^^). From this we get a modified relation 



s 



A tanh 



2T 



(14) 



with ?7 = Tin 



1 + exp ( f + In 



. The relation (HM captures the essential features of 



the purple curve in Fig. 2. 



B. A> Ec 

The green dashed line in Fig. 3 shows an example for the thermopower without parity 
effects for A > Eq- The arguments we have given in the preceding subsection for the 
thermopower at half-integer values of rix are valid also in this case. That is, Eq. ( !T3l) 
correctly describes the behavior of the function S{nx) also in this parameter range. Hence, 
this equation turns out to be the key for understanding the thermopower in NSN SET, i.e., 
a device with an electronic spectrum that is gapped around the Fermi energy. We further 
mention that for the thermopower without parity effects for the average energy close to 

9 



integer values of — n one has (^) ~ («„ + Un-i)/2 which yields 



b ~ — — , Tlx — n . (15) 



A=1.2£^, r=0.03£^^ 
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FIG. 3. Thermopower 5(n^) of the NSN SET with A = 1.2Ec and T = O.OS^c a) including 
all tunneling processes (purple solid line), b) without two-electron tunneling and without escape 
processes (green dashed line). While the green line is e-periodic the purple curve shows clear 
2e-periodicity and substantial qualitative changes with respect to the case without even-odd effect. 



If the gap exceeds the charging energy, two-electron tunneling starts to play a prominent 



role for the current- voltage characteristics of NSN SET in the parity regime T < T* |25| . 
Around odd-integer values of Ux there occurs a current peak which is due to a cycle of two- 
electron tunneling processes. Note that these processes do not have a gap. The thermopower 
is analogous to that of single-electron tunneling in NNN SET, the only difference is that we 
have to substitute the charging energy difference for two-electron tunneling, (^) ^ {un{nx) + 
u„_i(n^.))/2 which leads to a thermopower 5* — ('«« + '"n-i/(2eT). This result curiously 
coincides with the one in the absence of parity effects, Eq. (fT5|) although the dominating 
tunneling process is a different one. 

In the vicinity of the half-integer values of we observe the analogous effect to the case 
A < Ec'- The curve gets shifted towards the closest even-integer Ux without changing the 
slope. Scrutiny of the dominating current-carrying processes, e.g., for < 1/2 reveals that 
also here the largest rate is due to quasiparticle recombination on the island. However, if T < 
T* there are essentially no thermal quasiparticles. There is only a single unpaired electron 
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which is left from the pair-breaking tunnehng off the island. The quantitative description of 
this recombination process leads in full analogy to the conclusion of the previous subsection, 
namely that Eq. (HM describes this shift of the curve. The essence also here is that the 
current is due to tunneling of a single unpaired quasiparticle. 



V. FIGURE OF MERIT ZT 



Once we have studied the thermopower of the NSN SET it is an interesting question to 
investigate the thermoelectric efficiency of this device. This efficiency is quantified by 



ZT 



GvS^T GyS^T 



(16) 



where Ke and denote the electronic and lattice contribution to the heat conductance of 
the device. Our work focuses on a regime of extremely low temperatures for which is very 
small, hence we neglect it and k, = Kg. We compute k from the heat current Ig according to 
Eqs. ©-dH). 

An example of the results is shown in Fig. 4. For A ~ Ec we observe huge values of 
ZT which is rather uncommon, keeping in mind that typical values for materials reach the 
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FIG. 4. Figure of merit ZT for an NSN SET with A = 1.2Ec and T = OMEc a) including 
all tunneling processes (purple solid line), b) without two-electron tunneling and without escape 
processes (green dashed line). A strong enhancement of ZT is found in a range of Ux values where 
also the modulus of the thermopower reaches its maximum (cf. Fig. 3). 



order of 1, or for quantum dots the order of 10^. The question is why such high values are 
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possible for this system. A quick answer might be that, given that superconductors are bad 
heat conductors, it could be expected that by including the superconducting island ZT of a 
single-electron transistor would be enhanced. Let us try to give a more quantitative answer 
which takes into account our observations from the previous section. 

By transferring Matveev's idea for the interpretation of S as an average transport energy 
of the electrons according to Eq. ( !T2l) we can derive an expression that illuminates the 
meaning of ZT . In analogy with Eq. (IT^ we obtain K = {^'^)Gv / {e^T). This leads, 
together with Eq. (1161) and k = K — GyTS^, to the new relation 

- ■ 

Indeed, all additional factors cancel and ZT turns out to be the ratio of the squared average 
transport energy and the variance of that energy. This relation clearly indicates the strategy 
that needs to be used in order to increase ZT: The current-carrying electrons should be far 
from the Fermi energy while their energetic distribution should be as narrow as possible. This 



corroborates also the conclusion of Ref. 
for a high ZT value. 



19| that a 6 function in the spectrum is favorable 



Let us now apply Eq. (fT7j) to the NSN SET. We have already discussed that for the 
values where the strong enhancement of ZT is found, the current is carried by a single 
unpaired quasiparticle. The maximum energy of this quasiparticle is ~ A while its energy 
distribution is rather narrow: it is just given by the temperature T. Hence we expect 
ZT ~ A^/T^. In fact, this estimate has the correct order of magnitude. 



VI. DISCUSSION 

We have investigated the thermopower for a single-electron transistor with a supercon- 
ducting island for arbitrary ratios A/Ec and for temperatures both above and below the 
crossover temperature for parity effects T*. The results show the expected parity effects also 
in the functional dependence S{nx) of the thermopower on the gate charge n^, in particular 
2e periodicity. We have provided a discussion of the essential features of this functional 
dependence in terms of Eqs. (fT^ . (fT5|) . It is remarkable that the basis for this discussion is 
Eq. (fT3l) which was found already in Ref. (ll| . 

Apart from the thermopower we have also calculated the thermoelectric figure of merit 
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ZT. Unexpectedly we have found a strong enhancement of ZT compared to common values 
of this quantity. In order to understand our findings we have given a new interpretation 
of ZT in terms of Matveev's idea to represent thermoelectric quantities as moments of the 
energy distribution for the current-carrying electrons. It shows that large values of ZT can 
be obtained if the dominant transport mechanism occurs far from the Fermi level, and at 
the same time, has a narrow distribution in its energies. Clearly, the NSN SET is a system 
where these conditions can be achieved. 
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